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Galileon models arise in certain braneworld scenarios as modifications to General Relativity, and 
are also interesting as field theories in their own right. We show how the galileon model can be 
naturally generalized to include local gauge symmetries, by allowing for couplings to Yang-Mills 
fields. The resulting theories have at most second order spacetime derivatives in any representation 
of the gauge group, thereby avoiding Ostrogradski ghosts. We also extend the models to include 
curved space, and show how in that case we need to include non-minimal couplings between the 
galileons and the curvature tensors for the theory to retain its second order nature. 



I. INTRODUCTION 

Recently, there has been plenty of interest in galileon 
models, either as modifications of gravity or as effective 
field theories (see l2i] for detailed reviews). Galileons 
differ from conventional scalar fields in that they possess 
a Galilean field symmetry 

TT — > TT + 6^a;^ + c, fop, c = const., (1) 

hence the name This field symmetry was mo- 

tivated by the braneworld Dvali-Gabadadze-Porrati 
(DGP) model [4j], whose behavior within the Hubble hori- 
zon can be captured by a 4D boundary effective field 
theory where the General Relativity tensor modes are 
linearized and many distinct features of the model are 
carried by a galileon that encodes how the DGP brane 
is bent in the bulk Q . Remarkably, this generalized 
shift symmetry of the field is quite constraining in terms 
of constructing all possible Lagrangian terms, if we also 
require there are no higher than second order derivatives 
in the equation of motion. This second order restriction 
is important, as we require a well-defined Cauchy prob- 
lem and higher order derivative terms generally give rise 
to Ostrogradski ghosts (see, e.g. Q). In four dimensions, 
there are only 5 possible galileon terms Unlike the 
DGP model the general galileon model can give rise 

to a ghost- free self-accelerating branch Q . From the per- 
spective of effective field theory, the galileon model has a 
healthy energy gap between the classical non-linear scale 
and the quantum scale around a heavy background @ , 
analogous to that of General Relativity, and the couplings 
of the galileon interactions are not renormalized by loop 
corrections [1,0. 

As a surprising and interesting application, the galileon 
construction has proven crucial in building consistent 
massive gravity models. For a long time, a consistent 
massive gravity theory seemed unattainable, due to prob- 
lems such as the Boulware-Deser (BD) ghost Re- 
cently, inspired by the galileon construction, non-linear 
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massiv e g ravity without the BD ghost has been pro- 
posed [ll| - [l4| . where the helicity mode of the massive 
gravton is exactly a galileon. 

Due to these novel properties and appications, inter- 
est has spread to consider more general galileon models, 
motivated by their phenomenologies, 4D symmetries or 
braneworld scenarios. Since the galileon can be identified 
as a brane bending mode, multiple galileons may arise 
in a braneworld scenario with higher co-dimensions [l^ . 
General multi-galileon models and their phenomenolo- 
gies such as self-acceleration and self-tuning have been 
studied [l^. In [ll], p-form galileon models have been 
proposed where the galileon scalars are generalized to 
be differential forms. In certain braneworld scenarios, 
if a higher co-dimensional bulk has some global symme- 
try, or isometry, the resulting 4D effective galileon model 
may inherit this symmetry, leading to a galileon model 
with global internal symmetry 0, Il7l - [l9l | . Indeed, in this 
case, the galilean field symmetry and the internal sym- 
metry can both come from the bulk Poincare symmetry. 
The supersymmetric extension of the galileon model has 
been studied in (20| and an extension to the most general 
scalar field theory that is of second order in nature has 
been proposed in [2l| - |23j . 

The Galilean field symmetry ([T]) may originate from 
some underlying general covariance. This is indeed the 
case in the DGP model [1, Q. Therefore, the Galilean 
field symmetry may be ultimately relaxed, and we may 
consider a generalized symmetry, of which the Galilean 
field symmetry is some appropriate limit. An infrared 
completion of the galileon model has been considered 
along the lines of promoting the galilean field symmetry 
to that of a conformal symmetry Q, while Deffayet et 
al have constructed the covaria nt g alileon by directly co- 
variantising the galileon terms [24| . The general galileon 
model constructed by 4D symmetry considerations, as 
well as its conformal extension and direct covariantisa- 
tion, has been shown to arise in a probe brane setup 
where the brane and bulk are endowed with the Love- 
lock generalisation of General Relativity [1^. The gen- 
eral galileon model is the small field limit for a 3-brane 
in a 5D Minkowski bulk, while the conformal galileon 
arises for the case of an anti-de Sitter bulk. When the 
induced brane metric is allowed to fluctuate, the covari- 
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ant galileon is also recovered. More general brane and 
bulk settings have been considered, leading to more gen- 
eralized galileon field theories in 4D [2^ . 

On the other hand, the galileon models contain higher 
order derivative terms in the Lagrangian, which are sta- 
ble under quantum corrections [0, joj and may be used 
to construct solitonic solutions jlTl l27l|. In particular, an 
explicit soliton with a non-linear sigma model constraint 
has been obtained in the SO (4) symmetric galileon jTzj . 
Further discussions of the gravity and field theoretic as- 
pects of the galileon models may be found in (28l - l3l| . 
The galileon models have also sparked plenty of interest 
in cosmology, when considering both the late [s^l and 
early universe [ssj . 

As we have mentioned, a global geometric bulk sym- 
metry in certain braneworld scenarios can induce a global 
internal symmetry on the corresponding 4D effective 
galileon theory. However, generically, the bulk geometric 
symmetry may be realized locally, i.e., there is some un- 
broken general covariance alongside with possible isome- 
trics in the bulk. In this case, the resulting effective field 
theory would inherit a gauged internal symmetry. In- 
deed, from the formalism of a probe brane plus an in- 
duced metric of [25i] . this is what we would have when 
there is a bulk metric whose symmetries amount to re- 
placing the partial derivatives in the induced brane met- 
ric with gauge covariant derivatives. 

In this paper, we consider an important generalisa- 
tion of the galileon model, asking the question whether 
the galileons can be consistently coupled to gauge fields 
whilst maintaining the attractive features such as avoid- 
ance of Ostrogradski ghosts. In Section HTJ we briefiy 
review the galileon model with global internal symme- 
tries. We argue that the most important feature of the 
galileon model is that the spacetime indices of galileon's 
derivatives are appropriately anti-symmetrized and we 
shall keep this feature when gauging the galileon model. 
In Section [ml we construct the minimal gauge extension 
of the symmetric galileon model, focusing on the funda- 
mental representations of all the classical groups. We 
show that its equations of motion remain second order, 
naturally avoiding Ostrogradski ghosts. This is in con- 
trast to the gravitational covariantization of the galileon 
model, where Ostrogradski ghosts do arise unless extra 
terms with non-minimal couplings between the galileon 
and the curvature tensors are added to the theory. In Sec- 
tion |lVl we further promote the metric to be dynamical, 
in which case we show that non-minimal couplings be- 
tween the galileon and the curvature tensors are required 
in order to maintain the second order nature of the field 
equations. In Section |Vl we discuss non-minimal gauge 
covariantisations from the global symmetric galileon and 
gauged galileon models in other representations of the 
gauge group. We conclude in Section IVTl 



II. GALILEONS WITH GLOBAL SYMMETRIES 

First, we briefly review the galileon models with global 
internal symmetries. The galileon can be identified as 
the brane bending mode in certain braneworld scenarios. 
In a braneworld setup with several co-dimensions, the 
brane has several directions to bend, and the resulting 4D 
effective theory may contain several galileons. Neglecting 
the tadpole terms, the general A^-galileon model in 4D ^ 
can be written as [l^l 

5 

n=2 

(2) 

where ij,^^,...,^^ label the different galileons, Xi-^i2..A„ are 
free parameters of the theory and summations over the 
labeling indices zi,...i„ are assumed. Note that we de- 
fine the n = 2 case as Xi^i^d^^'^TTi^d^^'Ki^. All the down 
spacetime indices ^^v-jMn ^''^ anti-symmetrized. Since 
all the up indices are contracted with the down indices, 
the up indices can also be considered as being 

anti-symmetrized. Because of this anti-symmetrisation 
of the spacetime indices, we can freely "move" a par- 
tial derivative from one Tr^j. to another by integration 
by parts. This is a key feature that gives the galileon 
models their interesting properties such as Galilean field 
symmetry, absence of Ostragradski ghosts and the non- 
renormalisation theorem (e.g. 0, Another result of 
the anti-symmetrisation of the spacetime indices is that 
the indices of Xi^i2...in can be chosen as symmetric, i.e., 
= \iii2...i„)- By simple combinatorial counting, 
we can see that the number of free parameters of a gen- 
eral multi-galileon model increases very rapidly with N: 
j:L2iN + n-iy./n\{N-l)\^. 

However, if the bulk of the underlying braneworld 
scenario contains some global symmetry, the resulting 
galileon theory may inherit this symmetry and the al- 
lowed Lagrangian terms become more restricted [1, [13] . 
That is, the multiple galileons now form some represen- 
tation of a Lie group 

TT = (7ri,7r2,...,7rAr)^, (3) 

and, apart form the Galilean field symmetry, the La- 
grangian is invariant under a global symmetry transfor- 
mation 

TT-^e^'-^V, (4) 

where Ca's are constants across the spacetime, T"'s are 
the generators of the corresponding representation, and 
summation over the group indices is assumed. The fun- 
damental representation of a group is arguably the most 
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basic irreducible representation, which defines the corre- 
sponding matrix group. Now, if the galileon lives in the 
fundamental representation ^ of a classical group SU(N), 
SOfN) or Sp(N), the galileon Lagrangian is simply given 

C = -d^TT^d^'TT + A(9^7r^9[^7r)(a''9^7r1'9''5p]7r), (5) 

where ^ is the Hermitian conjugate and A is a free param- 
eter. That is, the theory has only one coupling constant. 
Mathematically, the Lagrangian is so simple because, in 
the fundamental representation of the classical groups, 
the only symmetric primitive invariant of the correspond- 
ing Lie algebra which the symmetric galileon lives in is 
S^j ) where the down and up group indices differen- 
tiate the fundamental representation and its conjugate 
respectively. Strictly speaking, for SO(N), there are also 
5^^ and Sij, but the fundamental representation of SO(N) 
is a real representation, so they can be treated as the 
same group tensors. For the fundamental representation 
of the exceptional groups, there may be extra Lagrangian 
terms. For example, there is another symmetric primi- 
tive invariant d*^'^ in the Eq Lie algebra [lH, implying 
that under this symmetry the Lagrangian will contain 
additionally a cubic term contracted with this d^^'^. 

More generally, the galileon may live in other repre- 
sentations of the internal group. The Lagrangian is more 
complicated for these cases. For example, in all pos- 
sible Lagrangian terms for the adjoint representations of 
SU(N) and SO(N) were presented. To have a practical 
taste, all possible Lagrangian terms for the adjoint rep- 
resentation of SO(N) with N > 5 are given by 

Tr($n$), Tr(<9$a[$a^$), 

Tr(a$a[$aa$9^$), Tr(a$a[$)Tr(59$a^*), 
Tr(a$a[$aa$9(9$a<^$), 

Tr(a$a[$)Tr(9a$aa$9(^$), (6) 

where we have suppressed the spacetime indices for sim- 
plicity and $ = ($/) is an iV X iV traceless hermi- 
tian matrix, which lives in the adjoint representation 
of SU(N), transforming under an SU(N) defining matrix 
U = iU/)NxN as $ ^ Um^i. 



III. GALILEONS WITH GAUGE SYMMETRIES 

In this section, we promote the internal symmetry of 
the symmetric galileon to be a local one. That is, we 
couple the galileon to a gauge field = ^J^T", where 
T"'s again are the generators of the symmetric galileon's 
representation. We will work in the general case where 



the gauge field is of the Yang-Mills type and the expres- 
sions will be written in terms of matrices that live in 
the representation of the corresponding Lie algebra. The 
Abelian theory can be obtained as a special case where 
is merely a real function of spacetime ^. 
We first establish our convention for the gauge field 
Af^ = A^T° and the gauge covariant derivative Z?^. We 
define the field strength as 

F^. = F;,T^ (7) 
= 2di^A,]+i[A^,A,] (8) 
= 2d[^Al^T^^r''-AlAlT^, (9) 

where /"'"^ are the structure constants of the gauge group 
satisfying 

[T'',T'']^ir^''T''. (10) 

We are looking for a Lagrangian that is invariant under 
a gauge transformation 

TT ^ n{x)TT ^ e*""(^)^°7r, (11) 
A^, ^A^Sl'^ + (12) 

where again tt is a column vector and Q;a(x)'s are func- 
tions of spacetime. A simple prescription for the gauge 
covariantisation is to promote the partial derivatives in 
the galileon theory to gauge covariant derivatives: 9^ — > 
Z?^. Later, we will also encounter gauge covariant deriva- 
tives in curved space. From a viewpoint of differential ge- 
ometry, we could simply call all these derivatives as the 
covariant derivative that takes into account all the fibres 
on a base manifold and changes its form according to the 
object it acts on. But for concreteness we choose to use 
different symbols for these different covariant derivatives. 

In Minkowski space, the action of the gauge covariant 
derivative on the galileon tt is given by 

= {d^ + zA^)7r = (9^ + iAp:'')Ti. (13) 

Note that the gauge coupling does not appear in the 
gauge covariant derivative but will appear in front of the 
gauge kinetic term. Since the field strength lives in the 
adjoint representation 

F^, ^ (14) 

the gauge covariant derivative on F^i^ is given by 

DpF^. = dpFp, + i[Ap, F^,]. (15) 

We will also use the left acting gauge covariant deriva- 
tives Dp — (Dp)^ , whose actons on tt^ and Ff^^, are 
defined by 

TT^Dp - nHdp - lAp) = dpn^ - mUp (16) 
Fp,Dp = DpFp,. (17) 



^ Of course, to construct a Lagrangian invariant for a complex 
representation, its conjugate also has to be included. 



^ An explicit example of a U(l) gauged galileon model can arise 
in 3D New Massive Gravity [3I . 
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Now, if the galileon lives in the fundamental represen- 
tation of a classical Lie group SU(N), SO(N) or Sp(N), a 
direct promotion from D/^i for the Lagrangian ([5]) 

leads to the following Lagrangian 

+ Xi ((7^'^7r)ti?[^7r) {iD''D,7r)^DPDp]Tr) 
^9g 

+ A4 (tt^D'^Di^tt'^ (n^DM'DPDp^Tr'^ , (18) 

where go is the gauge coupling constant and A4 may be 
called the galileon coupling constant. Note that in the 
Lagrangian (|18p both the up and down indices /i, p in 
the quartic tt term are anti-symmetrized. We will see that 
the equations of motion of this Lagrangian does not con- 
tain higher order derivatives. However, the Lagrangian 
psp , which we will refer to as the minimal gauging case, is 
not the only possible way to gauge the symmetric galileon 
model (O . We will discuss this later in Section |Vl Our 
conclusion remains the same if all the possible ways of 
gauging the Lagrangian ([5]) are considered. 

To facilitate the derivation of tt^'s and A'^'s equations 
of motion, we establish a few further relations regarding 
the gauge covariant derivative. The commutators of 

and with themselves are 

[D^,D,]TT^iFf,,n, n^[D^,D,]^nT^F^,. (19) 

Note that when acting on the field strength, the com- 
mutator of gauge covariant derivatives obeys a different 
rule: [D^, D^jFpa- = i[Fp^^, Fpa-]. To derive an equation 
of motion, we need to integrate by parts and use the 
Leibniz rule. In Minkowski space, this is with respect to 
the partial derivative and assumes integration over the 
total derivative 9^(...)^ vanish. In a gauge theory, the 
natural derivatives appearing in the action are the gauge 
covariant derivative. Nevertheless, as one may expect, it 
can be shown that the usual rules of derivatives apply for 
the gauge covariant derivative. Particularly, we have the 
following relations 

Jd\{7j^Dp)tl,'^{cl>^^) ^- JdWDM^{<l>^^)), (20) 

and 

(21) 

where 77, ip, </> and if take values on the same represen- 
tation as TT (a column vector in matrix form) but may 
carry gauge covariant derivatives and properly contracted 
spacetime indices. There are also similar relations where 
the up index is carried by 77, (j) or ip. As an aside, when 
r] and live in a representation different from that of (f> 
and if, these relations still hold, provided the gauge co- 
variant derivatives acting on these group vectors are de- 
fined appropriately in the corresponding representations. 



In short, we can perform integration by parts and use 
the Leibniz rule for gauge covariant derivatives in a way 
similar to that of the partial derivatives. 

With these rules established, we can easily derive the 
equation of motion for tt^ (tt's equation of motion is the 
hermitian conjugate of tt^ 's) . For a simpler presentation 
of the results, we define 

Jlt^^ - = Df^O^---7r, n-^^ =7rt...S^S^. (22) 

In matrix form, tt^'s equation of motion is given by 

£^'=£f+Xi£f, (23) 

where 

Sf - n'^^, (24) 

and 

+ n[/pn^n%] + n^i.n^n.^] . (25) 

There are terms with 3 or 4 gauge covariant deriva- 
tives acting on the galileon, but, thanks to the anti- 
symmetrisation of the spacetime indices, all the higher 
order derivatives exactly cancel each other. To see this 
clearly, we use Eqs. ([19]), which leads to 

- ^F'^PF^^pnn^u^^ + ^F-^^ni^n^n^] 

- \d^,{f%^)^'^fp^^^^ + ^n''[,(7rtF^'^)S,n^] 

- hv\p','^F^^F,p^'K + '-F^^pWWW^^ 

- '-Dy,{F%i,WW p^ + '-W^^^D'^iF,^^^) 

-^n^i^F.^'n,]. (26) 

Since F^y only contains first order derivatives, tt^'s equa- 
tion of motion does not have terms with higher than sec- 
ond order derivatives. 

Next, we also need to show that A'^'s equation of mo- 
tion does not have higher order derivatives. In deriving 
's equation of motion, the following corresponding re- 
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lations are useful, 
Lagran. variation 



EoM terms produced 



(27) 
(28) 

(29) 



EoM of A° 



(30) 



where 5a means variation with respect to A"^ , (j)f^ and 
ip^i, hve in the same representation as tt, and the RHS 
of the arrow is the equation of motion terms arisen from 
the Lagrangian variation of the LHS. These correspon- 
dences can be checked by exphcit calculations. Using 
these relations, as well as the commutator relation (IT9t . 
the equation of motion for A"^ in matrix form is given by 



S — Sp -\- S2 ^ A4^4 



(31) 



where 



£f = -^DPPp^ £f = -m^TT^ + zttH^ (32) 
9g 



and 



■i-i-i 

(33) 



Therefore, none of the equations of motion for the the- 
ory (|18l) has higher than second order derivatives. In 
other words, the galileon model can be naturally ex- 
tended to couple to gauge fields without introducing Os- 
trogradski ghosts. As we will discuss in Section |Vl the 
same conclusion holds for any representation of the gauge 
group in Minkowski space. 



IV. GAUGED GALILEONS IN A CURVED 
BACKGROUND 

We now allow the spacetime metric to fluctuate away 
from Minkowski space and construct a gauged galileon 
field theory in this new curved space background. This 
may be related to a braneworld scenario where the brane 
has non-trivial intrinsic geometry. 

At an arbitrary point of the curved spacetime, there 
is a tangent space as well as a Lie algebra of the gauge 



group. But they are different fibres on the base mani- 
fold, so the covariantisation of the galileon theory with 
gauge symmetries obeys the "comma-goes-to-semicolon" 
rule [13 : 

D^ipr. = {dp + iAp)iP^ VptP^ = {Vp + iAp)il;^, (34) 

where is the (tangent space) covariant derivative mul- 
tiplied by the identity matrix of the corresponding rep- 
resentation of the gauge Lie algebra. Now, the commu- 
tators of Vp and T>p = {Vp)'^ are given by 

[Vp,V,]^j^ = iFp,^^ + R%p,V, (35) 
V't-p^, V,] - i^^^Fp, ~ ^^''R%^., (36) 

where ip" lives in the same gauge representation as tt and 
R"^ Pfiu is the Riemann tensor of the curved spacetime. 
We have assumed the underlying gravity theory is torsion 
free. Note that the Riemann tensor term does not arise 
when acting the commutator on the spacetime scalars tt 
or TT^. 

In the fundamental representation of the classical 
groups, a direct use of the "comma-goes-to-semicolon" 
rule to the Lagrangian (fT8|) leads to the Lagrangian 



cs 



R ~ -LtriFp.F'^n - CDpTTyV^TT 
^9g 



+74 {{V^'n)^V^pn) {{V''V,7r)^VPVp]7:) 



(37) 



where R is the Ricci scalar and 74 is the galileon coupling 
in curved space. But this Lagrangian actually leads to a 
higher than second order equation of motion for tt^ , sim- 
ilar to the case of the covariant galileon studied in [23]. 
Specifically, there is one higher order derivative term ap- 
pearing in TT^'s equation of motion, namely, 

i(V[.i?p^'"')S|,|I^S^] = lv''G^.S>S^ (38) 

where we have used Bianchi's identities, Gpi, = Rpi, — 
Rgpy/2, and, analogous to Eq. ([22)) . we have defined 



V^-V,, • • • TT, 



TT^ • • • VV,. 



(39) 



So the Lagrangian (1571) is not the one we are looking for. 
We want a Lagrangian whose equations of motion con- 
tain at most second order derivatives and reduce to the 
Lagrangian (jl8p in the flat space approximation. This 
Lagrangian is given by 



C 



cs 



C 



cs 



-cs 

-C ' 



(40) 



where C'^ is a coupling between the Einstein tensor and 
the galileon 



L,Q — yy—gupu^ ^p^ ^ . 



74 



6 



(41) 



Now, the modified Lagrangian PO)) has a well-posed 
Cauchy problem, as we are about to check explicitly. 
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The equation of motion for tt^ is, in matrix form, 



and 



where 



and 



A.2 — ^ fj., 



(42) 



(43) 



t ^ 2 



p] 



p] 



- iG^.I?''(S,S^S'') ~ iG,,,I?^(S^S''Sp) 

7 1 ^ 

1 

- ^I?[.(^^%7r)vrtF''^S,] + ^E^^pin^ F^^.E,^ 

- Js^ip^t^p.^^^j^ + ^j.j^^spgps^^, 

+ - '-V[^{F%n)EPE%^ 

+ '-EypE^V^F,^^^) - '-E^pE'^F^^E,^ 

- l^'l.^^'R.Kf^x - \F'^PF^.p^E^E,y (44) 



Note that the index \ is not in the anti-symmetrisation. 
We can see that there are no higher order derivatives in 
TT^ 's equation of motion. 

We also have to check the equations of motion for A'^ 
and g"^ . The remarkable thing about the (counter) term 
(PT|) is that it also reduces all the higher derivatives in 
the equations of motion for A'^ and 5"'^. 

By explicit calculations, we can also establish a set of 
rules to facilitate the derivation of the equation of motion 
for A" , which are equivalent to those of (|27ll30p , with 
replaced by V^. Using these relations, for the modified 
Lagrangian (j40|) . A'^'s equation of motion is given by 



X"" = Xi 



X2 + 74<^4 



where 



Xp — —VFpa X2 
9g 



(45) 



(46) 



•^4 — '"[ct~ " p] + W<y[p~^ " p] 

1 
2 

7 <- 1 

-O r H ^''P^i- |7rt"P" , _ _o , "P Trt-Trtp^''^ , 



„ pl'PTJP'^ n R Aiyp7j/i77 

ig.[p7r(7rtF^)P^S^S,,] + ig,[,^S^^t^PPs^j 



ig^lpTTE" jyE'^E'' p] 



-G Tr~P~P~ 
g^Pcr"" " "P- 



(47) 



Again, the index \ is not in the anti-symmetrisation, and 
there are no higher order derivative terms. 

Finally, we show that the equation of motion for the 
metric g"^^ does not contain higher order derivatives. The 
full equation of motion for g°'^ is lengthy, and we do 
display it here. Instead, since we are only concerned 
whether there are higher order derivatives, it is suffice 
to only focus on the potentially dangerous terms, i.e., 
terms potentially having higher order derivatives, and 
check whether they indeed cancel each other. First, note 
that the conventional Lagrangian terms 



'9, 



(48) 



do not give rise to equation of motion terms with higher 
order derivatives. To derive g^^^s equation of motion 
for the remaining terms, note that the variations with 
respect to the metric g°'^ and ^/—g (rather than deriva- 
tives of g"^) do not give rise to higher order derivative 
terms. The only potential higher order derivatives come 
from the following variations 



(49) 



5Cl 



^-±^{5Rp, _ ]^g^P5R^pgp,)EPEpE^^E''. (50) 



where (JF^^ is the variation of the Levi-Civita symbol 
with respect to g°'^ (including g"^'s derivatives) and the 
down indices a and are not in the anti-symmetrisation. 
We will see that both of these two variations give rise 
to genuine higher order derivative terms, which exactly 
cancel each other. 

Neglecting the terms that obviously do not contain 
higher order derivatives, the equation of motion from 



7 



Eq. (gg]) includes 



1 



1 ^ 4- 



"/^" " ^P" 



2i!a[P^ "P"iy] 
" 2" "[A^^l/^pa p] 

-S^'^f^^i.^Q,]^/?. (51) 



The last 4 terms seem to have third order derivatives, 
but they are protected by the anti-symmetrisation of the 
spacetime indices. For each of these terms, two indices 
of the third order derivative are anti-symmctrized, thus 
we can use the commutator relations 



P5II36P to reduce 
them to lower derivative terms involving the gauge field 
strength or the curvature tensors. Therefore, for the vari- 
ation ()49p . the only non-reducible higher order deriva- 
tives are contained in 

1 ^„ ^ _„ 1 



■yg 



(52) 

On the other hand, neglecting the terms that do not 
contain apparent higher order derivatives, the equation 
of motion terms from the variation (1501) include 



12 "p," "p I- 



1 

12 
1 

12 

ff a ,3 ^ p 1/ 



■^^■^Q^''^/3pp 



1 

12 
1 

12 



1 fcp_^fcp„ 



"a^P" "P P 



a „~t^~'^" "P 



12 
1 

12 



12 



p^ 



^^^a^fjpii^'^ 



12 



1 

12' 



, „"P" p= ■=■" 

ap' — ' ^pp — ' 



12 
1 

12 
1 

12 
1 

12 



(53) 



Upon using the commutator 
fact that the indices „ and / 



12 
1 

12 
1 

' 12 
1 

12 

relations P5ll36p and the 
I are symmetric, each line 
of the last 5 lines can separately be reduced to lower or- 
der derivative terms involving the gauge field strength or 
the curvature tensors. Using the same technique and ne- 
glecting all the lower order derivative terms, the first 3 
lines and the second 3 lines can respectively be grouped 
together into the first and the second term of Eq. ([5^ 
but with the opposite sigh. Therefore all the higher or- 
der derivatives in g"^^s equation of motion exactly can- 
cel. There is no Ostrogradski ghost for the modified La- 
grangian 



Now, we want to compare the differences that arise 
between the gauge and gravitational covariantisations. 
When applied to the galileon model in both cases, the 
partial derivatives are replaced by the corresponding co- 
variant derivatives, which, unlike the partial derivatives, 
do not commute with themselves. Commuting two co- 
variant derivatives gives rise to a term involving either 
the gauge field strength (gauge case) or the Riemann 
tensor (gravity case). In the covariantized theories, the 
inherited galileon anti-symmetrisation of the spacetime 
indices leads to a significant reduction in the number of 
higher order derivatives in the equations of motion, and 
indeed eliminates all the higher order derivatives for the 
gauge case. However, This total elimination does not oc- 
cur in a gravitational covariantisation. The reason for 
the difference is due to the fact that, while the gauge 
field strength only contains first order derivatives, the 
Riemann tensor contains second order derivatives. 



V. NON-MINIMAL GAUGING AND GENERAL 
REPRESENTATIONS 

Now, let us discuss other ways to gauge the symmetric 
galileon with fundamental symmetries. We mentioned 
earlier that the Lagrangian JS]) is the only possible La- 
grangian for the galileon globally charged by the fun- 
damental representation of the classical groups. This is 
due to the fact that the partial derivatives commute with 
each other and the partial derivatives' spacetime indices 
are "dressed" with the anti-symmetrisation. Therefore 
the partial derivatives can be freely moved from one ni^, 
to another via integration by parts. However, this fea- 
ture is not kept intact in the gauge covariantized theories, 
both in flat and curved space. In these covariantized the- 
ories, we can still do covariant integration by parts and 
the anti-symmetrisation "dressing" still largely reduces 
higher order derivatives, but the covariant derivatives do 
not commute and their commutators give rise to terms 
with the galileon coupled to the gauge field strength or 
the curvature tensors. 

Indeed, in the fiat space case, the following terms 

{{D^7r)^D[^D''7r) {{D^ii)^ DP Dp^n) + h.c, (54) 
{{D^7r)^D[^D'^7r) {{D^DPtt)^ Dp^n) , (55) 

where h.c. means the Hermitian conjugate of the previous 
terms, are different from the A4 term of (|18p . Neverthe- 
less, thanks to the anti-symmetrisation "dressing" , these 
terms do not give rise to higher order derivative terms in 
the equations of motion, as one may explicitly check. In- 
deed, by integration by parts, these terms can be reduced 
to the A4 term of (1181) plus a number of terms involving 
non- minimal couplings to the field strength f . The sit- 
uation for the curved space case is similar. We will have 
the curved space counterpart of the terms (IM|) and ([55)) 
(with Z?p replaced by I?p). Apart from the non-minimal 
couplings to the gauge field strength, non-minimal cou- 
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plings to the curvature tensors will also arise from reduc- 
ing the curved space counterpart of the terms and 
([55| to the 74 term of ([40l) . 

For representations other than the fundamental of the 
classical groups SU{N), SO{N) and Sp{N), as we men- 
tioned, more complicated galileon terms, such as those 
of (l6|), are allowed in the Lagrangian with global inter- 
nal symmetry. When constructing the gauged theory, 
we will also have the choice of minimal gauging or non- 
minimal gauging. Nevertheless, in flat space, the gauged 
galileon theory, minimally or non-minimally gauged, do 
not have Ostrogradski ghosts in any representation of 
any gauge group, including the exceptional Lie groups. 
As one may have realized from our explicit calculations 
with the fundamental representation, the reason for this 
is that: In the equations of motion of both the galileon 
and the gauge field, the gauge covariant derivatives are 
"dressed" by the anti-symmetrisation of the spacetime 
indices and the most apparent higher order derivatives 
are fourth order; After using the commutator relations of 
the gauge covariant derivatives, there are only first order 
derivatives of the gauge field strength, which itself only 
has first order derivatives; Therefore in the equations of 
motion the most one can get is second order derivatives. 

VI. SUMMARY AND OUTLOOK 

Although originally proposed as modifications to Gen- 
eral Relativity on very large scales, the galileon models 
have taken on their own lives as effective field theories, 
due to their interesting field-theoretic properties. Var- 
ious generalisations of the original galileon model have 
been proposed recently. In particular, galileon models 
with global internal symmetry have been presented and 
explicitly investigated in a probe brane setup with bulk 
isometrics. In this paper, we have constructed for the 
first time the galileon models where the internal sym- 
metry becomes local, i.e., the galileon is coupled to a 
Yang-Mills gauge field. This gauge covariantisation will 
break the Galilean field symmetry, akin to the gravita- 
tional covariantisation of [23|. But in the limit where 
the local symmetry reduces to a global one, the galileon 
model with a global symmetry will be recovered. We have 
shown that this generalisation is natural in the sense that 
the Ostrogradski ghosts remain absent after the gauge 
covariantisation, despite apparent higher order deriva- 
tives in the Lagrangian. Note that it is indeed non-trivial 
for an apparent higher derivative Lagrangian term to be 
free of Ostrogradski ghosts. Although we only explic- 
itly focused on the minimal gauging of the classical Lie 
groups' fundamental representation, we have shown that 
the absence of Ostrogradski ghosts is generally true for a 
gauged galileon theory in fiat space in any representation 
of the gauge group. This is essentially because the gauge 
field strength only contains first order derivatives, while 
the Riemann tensor contains second order derivatives. 



We also couple the gauged galileon to gravity. For the 
naive covariantization, higher order derivatives do arise 
in the equations of motion. However, by adding an extra 
non-minimal coupling involving the Einstein tensor, the 
revised Lagrangian has a well-defined Cauchy problem, 
analogous to that of the covariant galileon 24] . 

We have argued that the gauged galileon can arise 
from a braneworld scenario, particularly from a probe 
brane setup. Galileon models with global internal sym- 
metry arise when the bulk has maximal isometrics, such 
as a Minkowski bulk with the induced brane metric 
V = ^X^/^x^'^XB/^x■'GAB = T^M-^ + a^Tr'a^TT,, where 
X^{x) is the embedding function and Gab is the bulk 
metric 0]. However, generally, the bulk may have some 
unbroken general covariance, as well as possible isome- 
trics. In this case, the brane induced metric may obtain 
a form like /i^i^ = g^j^ -|- {D^tt'^)'^ D^Tii, which leads to 
the gauged galileon in the fundamental representation 
we have explicitly considered. Simple counting of the de- 
grees of freedom indicates that this kind of ansatz is gen- 
erally obtainable, as there are more bulk metric compo- 
nents than the embedding equations needed to be satis- 
fied, even for the Abelian case. It is interesting to identify 
the underlying symmetries of this braneworld setup and 
construct the corresponding bulk metric, which we leave 
for future work (Recently, Goon et al have explicitely 
realised our conjecture [37|-)- 

We know gauge fields play a vital role in many self- 
consistent field theories and generally arise in stingy sce- 
narios with extra dimensions, so galileon models must 
also be able to accommodate internal gauge symmetries. 
The goal of this paper is to demonstrate how this ac- 
commodation can be successfully realized. In we 
explicitly constructed a solitonic solution in the global 
SO (4) galileon model. The obtained gauged gahleon 
terms may add extra features in constructing inter- 
esting non-perturbative field theory solutions such as 
monopoles, textures and cosmic strings. The presence 
of such Abelian and non- Abelian gauge fields can also 
have dramatic consequences in the associated cosmology, 
ranging from the generation of primordial or cosmic scale 
magnetic fields to the formation and evolution of cosmic 
strings. We may also consider the gauged galileon mod- 
els in curved space as modified gravity models, in which 
case the gauge field provide an extra vector mode, so we 
will have a scalar- vector-tensor theory. The consistencies 
of these models in different representations of the gauge 
group are yet to be investigated. 
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